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Abstract 

We study a model for polaron-like charge transport mechanism along DNA molecules 
with emphasis on the impact of parametrical and structural disorder. Our model 
Hamiltonian takes into account the coupling of the charge carrier to two different 
kind of modes representing fluctuating twist motions of the base pairs and H-bond 
distortions within the double helix structure of A— DNA. Localized stationary states 
are constructed with the help of a nonlinear map approach for a periodic double 
helix and in the presence of intrinsic static parametrical and/or structural disorder 
reflecting the impact of ambient solvent coordinates. It is demonstrated that charge 
transport is mediated by moving polarons respectively breather compounds carry- 
ing not only the charge but causing also local temporal deformations of the helix 
structure through the traveling torsion and bond breather components illustrating 
the interplay of structure and function in biomolecules. 
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1 Introduction 



Electronic transport through DNA is crucial for its biological functions such as 
the repair mechanism after radiation damage and biosynthesis. As proposed in 
[1] electron transport through DNA proceeds along a one-dimensional pathway 
constituted by the overlap between 7r— orbitals in neighboring base pairs. 
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With the perspective to possible applications in molecular electronics bioma- 
terials are supposed to play an important role because they may offer the ulti- 
mate way to miniaturization and their attributed electrical transport proper- 
ties are expected to be markedly different from those of traditional macroscopic 
conductors [2]. In this context enormous effort has been directed towards in- 
vestigations of conductivity of DNA viewed as a promising conduit for charge 
transport on nanoscales [3]. So far the findings concerning conductivity of 
DNA have been controversial. There are the results of recent measurements 
indicating that DNA behaves as a well conducting one-dimensional molecular 
wire [4]- [10]. In contrast it was reported that DNA is insulating [11] and for 
short oligomers built up from base pairs of the same type semiconductivity 
was observed [12]. This qualitative discrepancy is shared also by theoretical 
findings so that it remains unclear whether DNA is conductive or not. For ex- 
ample a theoretical investigation based on first principle calculations suggests 
the absence of dc-conductivity in DNA [13]. On the other hand there stand 
the studies using various assumptions for modeling the DNA structure which 
focus on different aspects such as the influence of aperiodicity, temperature 
driven fluctuations and aggregation effects on the energetic control of charge 
migration which have all come to the conclusion that DNA is a conductor 
[14]-[18]. 

The detailed understanding of the conduction mechanism of DNA remains still 
an unsolved task. Different attempts to model the charge transport of DNA 
were based on transport via coherent tunneling [1] , classical diffusion under the 
conditions of temperature-driven fluctuations [16], incoherent phonon-assisted 
hopping [19], [20], variable range hopping between localized states [17] and 
charge carriages mediated by polarons [21], [22] and solitons [23]. 

In the present work we use for the theoretical description of charge transport 
in DNA a nonlinear approach facilitating the concept of polaron respectively 
breather solutions. The structure of the bent double helix of A— DNA is mod- 
eled by a sterical network of oscillators in the frame of the base pair picture 
[24], [25], [26] taking into account deformations of the hydrogen bonds within a 
base pair and twist motions between adjacent base pairs. The electron motion 
is described by a tight-binding system. It is assumed that the electron motion 
is predominantly influenced by vibrational modes of the double helix. The 
nonlinear interaction between the electron and the vibrational modes causes 
the formation of polarons respectively electron-vibron breathers supporting 
charge localization and transport. 

The aim of the current investigation is to elucidate the role played by static 
parametrical disorder and complex structural effects due to irregular devi- 
ations from the ordered helical shape of the double helix (reflected in ran- 
dom arrangements of the positions of the bases) with regard to the charge 
transport in DNA. There arises at least one interesting question, namely: Is 
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long-lived and stable charge transport in DNA achievable at all under the 
impact of disorder ? To tackle this problem in our paper we exploit analytical 
and numerical methods from the theory of (nonlinear) dynamical systems in 
combination with the theory of (linear) disordered systems. The paper is or- 
ganized as follows: In the first section we describe the model for the charge 
transport along the steric structure of the bent DNA double helix. The sec- 
ond section is concerned with the construction of polaron respectively static 
electron-vibron breather solutions of the nonlinear lattice system when para- 
metrical and structural disorder effects are taken into account. In section 4 
we discuss the mobility properties of the polarons respectively electron-vibron 
breathers and their transport efficiency in dependence on the degree of disor- 
der contained in the on-site energies and/or in structural arrangements of the 
helix. In a subsequent section 

we invoke the Floquet analysis to explain also the mobility of the polarons and 
breathers in terms of their spectral properties. Finally, we present a summary 
and discussion of the results. 



2 Model Hamiltonian for charge transfer along DNA 

DNA chains are quite flexible and exhibit structural fluctuations influencing 
the transport properties. Our model is designed to cover the basic features of 
the DNA double helix structure needed for a proper description of the charge 
transport dynamics under the impact of fluctuating base motions. Despite the 
fact that the charge transport proceeds across a one-dimensional channel of 
base stacks spanning along the sugar-phosphate backbone so that DNA can 
effectively be viewed as a one-dimensional molecular wire [27] some signifi- 
cant details of the three-dimensional DNA structure have to be incorporated 
in the model. For the physical description of DNA molecules their detailed 
chemical structure can be neglected. They can be simply considered as bent 
double-stranded systems with an adequate number of parameters sufficient 
for modeling the considered physical process. In the frame of simple DNA 
models such as the base pair picture the bases have been treated as single 
nondeformable entities. The helicoidal structure of DNA is then conveniently 
described in a cylindrical reference system where each base pair possesses two 
degrees of freedom, namely a radial variable measuring the transversal dis- 
placements of the base pair (that is deformations of the H-bond) and its angle 
with a reference axis in a plane perpendicular to the helix backbone which de- 
fines the twist of the helix [26] . In Fig. 1 we depict schematically the structure 
of the DNA model. 

With relevance to charge migration the most influential fluctuating motions 
of DNA originate from tranverse vibrations of the bases relative to each other 



3 



Fig. 1. Sketch of the helicoidal structure of the DNA model, the bases being rep- 
resented by bullets. Geometrical parameters Ro, 6q, Iq and the radial and angular 
variables r n and 9 nn -i, respectively are indicated. 

within each base plane, i.e. the base pair distance stretchings/compressions 
and the torsional variations of the helicoidal twist dominating longitudinal 
(acoustic strand) phonons [16] because the sugar-phosphate chains are rather 
stretch proof [27]. Therefore it is justified to discard the impact of longitudinal 
vibrations and the considered motion is restricted to the base pair planes 
[15]. Furthermore, the angular twist and radial vibrational motions evolve 
independently on two different time scales and can be regarded as decoupled 
degrees of freedom in the harmonic treatment of the normal mode vibrations 
of DNA [28]. 

The Hamiltonian for the electron migration along a strand in DNA consists 
of three parts 

H = H e i + H ra d + H twist , (1) 

with H e i representing the part corresponding to the particle charge transport 
over the base pairs, H ra d describes the dynamics of the H-bond vibrations and 
H twist models the dynamics of the relative twist angle between two consecutive 
base pairs. The electronic part is given by a tight-binding system 

H el = J2 E n \ C n\ 2 ~ K.n-1 (<C„_i + C n C* n _ x ) . (2) 
n 

The index n denotes the site of the n— th base on one of the two strands and 
c n determines the probability that the electron occupies this site. E n is the 
local electronic energy and V nn _i stands for the transfer matrix element which 
is responsible for the transport of the electron along the stacked base pairs. 
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The vibronic part H raa - takes into account radial displacements of the base 
units from their equilibrium positions within the base pair plane connected 
with deformations of the hydrogen bonds linking two bases. The corresponding 
part of the Hamiltonian treating the bond vibrations classically and harmon- 
ically is given by 



H ra d — X! 



(K) 2 M n tf r rl 
2M, 



n 



(3) 



with the momenta (Pn) 2 conjugate to the radial coordinates r n measuring 
the deviation from the equilibrium length of a hydrogen bond along the line 
bridging two bases of a base pair, Q r is the frequency and M n is the reduced 
mass of a base pair. 



The twist motion part reads as 



H tmst = £ [ + J ""«5°" 1 1 ■ ( 4 ) 



2 J, 



n 



where 9 nn -i is the relative angle between two adjacent base pairs quantifying 
displacements from its equilibrium twist angle 9 and J n is the reduced moment 
of inertia. 



We emphasize that in the context of charge transport in DNA the nonlin- 
ear dynamics of large amplitude structural transitions are not of interest be- 
cause charge transport processes are thought to be mainly facilitated by soft 
ffuctuational vibrational modes allowing the harmonic approximation of their 
dynamics. 

The interaction between the electronic variable and the structure variables r n 
and nn -i stems from the parameter dependence of the electronic parameters 
E n and V nn -i. The diagonal term expressing the most efficient coupling is of 
the form 

E n = E° n + kr n , (5) 

and takes into account the modulation of the on-site electronic energy i?° 
by the radial vibrations of the base pairs. In turn the actual charge occupa- 
tion has its impact on the local radial distortion of the helix. As the transfer 
matrix elements V n n-\ are concerned we assume that they depend on the 
three-dimensional distance between two consecutive bases along a strand in 
the following fashion 

K.n-1 = V (1 - a4n-l) ■ (6) 

The quantity a regulates how strong V nn -\ is influenced by the distance and 
the latter is determined by 
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dnn-i = {a 2 + (R + r n f + (R + r n _ x ) 2 

- 2(^0 + r n ){R + r n -i) cos(# + 6 nn ^)} 1/2 - l , (7) 

with 

l = yfa 2 + U% sin 2 (^ /2). (8) 

Expanding the expression (7) up to first order around the equilibrium positions 
gives 

R 

dnn-i [( 1 - cos0 o ) (r n + r n -i) + sin6> R 9 nn -i } ■ (9) 

'0 

Realistic parameters for DNA molecules are given by [26], [27]: a = 3AA, 
R « 10A, fl = 36°, J = 4.982 x 10~ 45 kg m 2 , tt e = [0.526 - 0.744] x 10 12 s _1 , 
fi r = 6.252 x 10 12 s _1 , V ~ 0.1 eV and M = 4.982 x 10- 25 %. 

We scale the time according to t — > t and introduce the dimensionless 
quantities: 



M 2 - r K ~ El 

~r n , k n = , , £° = ^ (10) 



V n 5 ,v n / — j -^n t r 



Subsequently we omit the tildes. 

With the use of the expression (9) the equations of motion derived from the 
Hamiltonian (2)-(4) read as 



irc n = (E n + kr n )c n 

- (1 - a d n+ l,n) C n+ i - (1 - a dnn-i) C n _i (12) 

f n = ~r n - k \c n \ 2 - a^(l - cos 9 ) 



x 



{ K+lCn + C n+1 C* n ] + [<C„_i + CnC^J } (13) 

R 2 

9 n n-l = -n 2 9 nn -l ~ OiV-± Sin^^c,,.! + CnC*^], (14) 

'0 

and the quantity r = hQ r /V appearing in Eq. (12) determines the time scale 
separation between the fast electron motion and the slow bond vibrations. 
We remark that in the limit case of a = and constant E n = E the set 
of coupled equations represents the Holstein system widely used in studies of 
polaron dynamics in one-dimensional lattices [29]. Furthermore in the linear 
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limit case emerging for a = k = the Anderson model is obtained for random 
£°[30],[31]. 

Regarding the dynamical energy exchange between the various degrees of free- 
dom we state that due to the large differences between the electronic and vi- 
bronic frequencies we are able to prove rigorously that the energy exchange is 
suppressed on time scales growing exponentially with the ratio of the frequen- 
cies. Corresponding bounds on the energy exchange suppression are derived 
with the help of Nekhoroshev-like arguments [32], [33]. Furthermore applying 
the concept of the anti-integrable limit [34] we are able to prove the existence 
of breather solutions for the coupled system (12)-(14). The mathematical de- 
tails are delegated to a forthcoming article [35]. 

The values of the scaled parameters are given by r = 0.2589, f2 2 = [0.709 — 
1.417] x 10~ 2 , V = 0.0823, R = 34.862 and l = 24.590. In our model study 
we treat the electron-mode coupling strengths k and a for which no reliable 
data are available as adjustable parameters. In the forthcoming investigations 
we adopt the values of these free parameters such that not too strong defor- 
mations of the helix result which is essential for our used harmonic treatment 
of the dynamics of the structural coordinates. 



3 Localized polaron-like states 



Being interested in the simulation of a nonlinear charge transport mechanism 
in DNA we embark as a first step on the construction of localized station- 
ary solutions of the coupled system (12)-(14). Exploiting the fact that the 
adiabaticity parameter r is small, i.e. the fast charge transport along the 
base stacks and the slow bond vibrations and the even slower rotational twist 
motions evolve on distinct time scales, the inertia in Eqs. (13) and (14) are 
negligible. Therefore, one can solve in the adiabatic limit the resulting static 
equations which results in the instantaneous displacements 



r n = -k \c n \ 2 - a ^ (1 - cos6» ) 

X { K+l C n + C„+l<] + [<C n _i + C n c;_J } , (15) 

ol V -R 2 

0nn-l = p^ - ~T~ Sin6 [<C n _i + C n C* n _^ . (16) 

Upon insertion of (15) and (16) into Eq. (12) one obtains a nonlinear discrete 
Schrodinger equation for the electronic amplitude 
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lTC n = 



El 



k \c n \ 



ka ^ (1 — cos^o) 



X { [ C n+l C n + c n+lC„] + [c n C n _i + C n C„_ 1 ] j 
— ( 1 + b n+ \ n ) C n+ i — ( 1 + b nn -i ) C n _i , 



(17) 



with 




+ a (1 - cos6» ) ( [< +1 c„ + c„+ic*] 

+ 2 [ C n C "-l + C n C n-l] + [c*_ 1 C n _ 2 + C„-lC*_ 2 ] ) 

+ 7^2~ ^0 Sin2 °0 [ C *n C n-l + ]} (18) 



In order to search for stationary localized solutions (polaron states) we sub- 
stitute c n = $ n exp[— % Et/r] in Eq. (17) and obtain the difference system 



k a — — (1 — cos Oq) 



— ( 1 + B n+ \ n )$„ +1 - (1 + B n „_i )$„_!, 



(19) 



where -B nn _i is given by 



Snn-1=« (^) < U 



COS O ) 



^[|$n| 2 +|<fn-l| 2 ] 
-n-0 



+ « (1 - COS O ) ( [$n+l$n + $ B +1^] 

+ ^ i? 2 sin 2 [$;$„-i + ^-j} 



(20) 



The ground state of the system (19) is computed with the numerical map 
method outlined in [36], [37]. 

In Fig. 2 we depict typical patterns of polaron states. First we discuss the 
ordered (periodic) case arising e.g. for synthetically produced DNA molecules 
consisting of a single type of base pairs (e.g. poly(G)-poly(C) DNA polymers) 
surrounded by vacuum which provides realistic ideal conditions for the consid- 
eration of DNA conductivity in experimental and theoretical studies [12]. The 
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Fig. 2. Profiles of the stationary polaron states for the ordered (periodic) case. 

(a) The electronic amplitudes \& n \ 2 - 

(b) The static radial displacements r n . 

(c) The static angular twists 9 n . 

electronic contribution is shown in Fig. 2 (a). In the ordered case the electronic 
wave function is localized at the central lattice site and the amplitudes decay 
monotonically and exponentially with growing distance from this central site. 
We remark that the spatial extension of the polaron decreases with increased 
coupling strength k and there is a smooth transition from large to small po- 
larons. Analogously, enlarging the off-diagonal coupling strength a leads also 
to enhanced degree of localization. The size of the polarons plays a crucial role 
for their mobility in the sense that for large up to medium polarons coherent 
motion can be activated whereas small polarons are immobile due to their 
pinning to the discrete lattice [38], [39]. The associated patterns of the static 
radial and angular displacements of the bases are shown in Fig. 2 (b) and (c), 
respectively. Like the electronic wave function the radial and angular displace- 
ments are exponentially localized at the central lattice site (base pair). Notice 
that due to the overall minus sign of the radial and angular excitation patterns 
the H-bridges get compressed while the helix experiences a local unwinding 
around the occupation peak of the localized electron. 

For a more realistic study we take static diagonal disorder in the on-site elec- 
tronic energy E® into account. The random values caused for example by the 
inhomogeneous broadening of the sites of distinct ion pairs with different ener- 
gies were simulated by random potentials [-AE, AE] with mean value E = 
and different mean standard deviations AE. (Note that any Eqc u term on the 
r.h.s. of Eq. (12) with constant Eq can be eliminated by a gauge transformation 
c n -> exp(-iE t/r)c n .) 
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Fig. 3. Breather motion along the DNA in the ordered case. Motion is initiated by 
an initial kick of the velocities {r n } directed along the corresponding pinning mode 
component. The kick amplitude is A r = 0.25. 

(a) Uniformly moving electronic breather. 

(b) The associated moving breather in the radial displacements. 

(c) The angular breather contribution consisting of a small amplitude breather 
pinned at the starting site and a relatively large amplitude moving breather com- 
ponent propagating in alliance with the electron and radial breathers. 
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Generally the localized excitation patterns do not change qualitatively when 
diagonal disorder is taken into account. Compared with the ordered case one 
merely observes a shift of the excitation peak away from the central lattice 
site depending on the respective realization of the random potential. However, 
since in the presence of diagonal disorder the translational invariance of the 
lattice system is broken the localized amplitude pattern is no longer reflection 
symmetric with respect to a distinguished lattice site as opposed to the corre- 
sponding symmetric localized state of the ordered system. Moreover we find 
that the combined effect of the two localization mechanisms, viz. nonlinear 
polaron and linear Anderson-mode formation respectively, leads to enhance- 
ment of the degree of localization compared to the case when only one of the 
two mechanisms acts. 

Furthermore, real DNA molecules exhibit random structural imperfections of 
their double helix caused e.g. by the random base sequence, the deforming im- 
pact of the chemical surroundings when DNA gets buffeted by water molecules 
and/or the varying hydrophobic potential of the base pair interactions de- 
pending on the ambient aqueous solvent that may leave the involved helix 
structures in irregularly distorted shapes. Accordingly, the structural disorder 
is incorporated in our model by randomly distributed equilibrium positions of 
the bases so that the resulting irregular helical DNA matrix deviates from the 
perfectly regular helix structure. 

To be precise, we consider randomly distributed torsional coordinates Q nn -\ 
and radii r n . Since the positions (actually, the distances) govern the value of 
the corresponding transfer matrix elements V n n-\ — (1 ~~ arf n u-i) nondi- 
agonal disorder is induced in the latter at the same time. Random transfer 
integrals along the helix were used to represent random base pair sequences 
for which the charge transport along DNA was discussed recently in a model 
relying on localized electron states [15]. 



4 Charge transport mediated by mobile polarons and breathers 

In this section we study the charge transport supported by mobile polarons 
respectively breathers propagating along DNA. To activate polaron motion 
an established method is provided by proper stimulation of localized internal 
normal modes of certain shape and frequency which are also called pinning 
modes [40]. The identification of possibly existing pinning modes is related to 
the linear stability analysis of the polarons. In order to investigate the stability 
of the polarons with respect to linear perturbations we derive the linearized 
system for which we impose small perturbations on the stationary polaron 
solution (<&„, r^\9^) in the form: 
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c n (t) = [$ n + Sc n (t)] exp(-iEt/r) 
r n (t) = rW + 6r n (t) 

e n (t) = ew + 8e n (t), 



(21) 
(22) 
(23) 



where 8c n , 5r n and S8 n are small time-dependent variables. Substitution of 
(21)-(23) into Eqs. (12)-(14) and retaining only the terms linear in the per- 
turbations results in the system of tangent equations. The stability analysis 
is then performed using Floquet theory. 

To this end we integrate the linear system of the tangent equations over one 
period T = 2nr/E yielding the Floquet map 



with the Floquet matrix F determining the evolution of the initial set of vari- 
ables x(0) = (Re5c n (0),Im5c n (0),5r n (0),5f n (0),56 n (0),59 n (0)) after times 
nT. To ensure linear stability all eigenvalues exp(i ■?/>„) of the Floquet ma- 
trix F have to lie on the unit circle, i.e. the Floquet arguments ip n have to be 
real- valued. Intensive numerical investigations in a wide parameter range have 
proved that all the polarons gained from the map method are linearly stable. 

Moreover the Floquet analysis supplies the frequencies of the normal modes 
via the relation uj^ = u[ip/(27i) + m] with ip expressed in rad and m is an 
arbitrary integer number. We identified the pinning modes for given sets of 
parameters. According to [40] we initiate the motion of polarons through suit- 
able perturbations of the velocity variables r n (t) and/or 6 n targeted in the 
direction of the pinning mode, that is we use for the numerical integration of 
the system (12)-(14) the following initial conditions 



with the normalized momentum parts £ r> of the pinning mode and \ r fi are 
the perturbation strengths (for details see [40]). 

Ordered case 

In Fig. 3 we present results for the activation of breather motion along DNA 
in the ordered case. We integrated the set of coupled nonlinear equations 
(12)-(14) with a fourth-order Runge-Kutta method and the accuracy of the 
computation was checked through monitoring the conservation of the total 
energy as well as the norm J2 n \ c n(t)\ 2 = 1- The applied kicking strength 
is A r = 0.02 corresponding to an amount of radial kinetic energy of the or- 
der of 200 meV. Interestingly, we found that there exists no angular pinning 
mode component (see further below) so that alone the excitation of the 
radial component £ r produces long-lived stable polaron and breather motions. 
Hence, our findings suggest that higher frequency radial fluctuational modes 



Sx(nT) = F n 5x(0) , 



(24) 



{ $ n , 0, r<°\ 0, 0<°>, } + {0, 0, 0, A r £ r , 0, A<&} , 



(25) 
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of DNA instigate the propagation of localized structures along DNA rather 
than low frequency twist fluctuational modes. The Fig. 3 (a) demonstrates 
that the electronic component of the polaron moves with uniform velocity 
along the lattice maintaining its localized shape throughout the journey. In 
the same manner the radial polaron part travels as a localized wave along the 
strand so that the electron propagation is accompanied by a local and tem- 
poral contraction of the radii. However, the angular static polaron component 
splits up into two parts when the radial pinning mode is initially imposed. At 
the starting position there remains a standing breather of comparatively large 
amplitudes which reverse periodically their sign in the course of time. Thus 
the helix experiences alternately a local winding and unwinding around its 
central base pair. In addition a second breather of relatively small and exclu- 
sively negative amplitude propagates in unison with the electronic and radial 
polaron components. The result is that the region of the helix which is tra- 
versed by the electron gets locally untwisted. As the quantitative assessment 
of the charge transport is concerned we note that the breather travels across 
seventy sites (base pairs) in a thousand (dimensionless) time units meaning 
that physically the charge conduction across the base pairs of DNA over a 
distance of 238 A takes ~ 0.16ns leading to a higher transport rate than a 
typical one found for protein conductivity [15]. This underlines the status 
ascribed to DNA molecules as promising candidates establishing fast and ef- 
ficient charge transport in (albeit only theorized yet) molecular electronical 
engineering utilizing them as molecular wires. 

Static diagonal disorder 



In the presence of diagonal static disorder, that is when the on-site electronic 
energies E® are represented by random numbers distributed in the interval 
\E®\ < AE with width AE, we find that mobile breathers exist up to a 
critical degree of disorder AE crit and conductivity persists. (Let us stress that 
it was indeed experimentally found that the random sequence A— DNA is an 
electrical conductor [9].) Beyond this critical value continuous propagation of 
the breathers along the lattice is prevented. 

This behavior is illustrated in Fig. 4 displaying the spatio-temporal evolution 
of the electronic amplitude |c n (t)| 2 . While for a relatively small amount of 
disorder AE = 0.05 the electronic breather (together with its radial and tor- 
sional counterparts) is able to travel uniformly along the lattice for a ten times 
increased AE = 0.5 the breather motion is restricted to a neighborhood of the 
initially excited position around which it wanders in an oscillatory manner. 

Regarding mobility there act two competing mechanisms in the combined 
disordered nonlinear system, namely a linear Anderson-type one related with 
the effects of disorder which attempts to pin the modes (bear in mind that 
the linear Anderson-modes are supposed to be immobile) and opposed to this 
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Fig. 4. Electron breather motion along DNA with diagonal disorder for two different 
values of the width AE of the random distribution of the on-site electronic energy 

(a) Width AE = 0.05. The electron breather moves with uniform velocity along a 
strand. 

(b) Width AE = 0.5. Confined electron breather with positions oscillating itiner- 
antly around the initial lattice position. 

there operates the nonlinear mechanism utilizing the pinning mode to activate 
motion. The relative strength of these two mechanisms AE/X r determines then 
breather mobility. In order to summarize the results for the breather mobility 
in the presence of diagonal disorder we display in 5 the temporal behavior of 
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the first momentum of the electronic occupation probability defined as 

n(t)=Y,n\c n (t)\ 2 , (26) 



30 - 

20 - 




Fig. 5. Diagonal disorder. The position of the center of the electron breather as 
a function of time and for different amounts of disorder. Full line (Ordered case), 
short dashed line (AE = 0.025), dashed-dotted line (AE = 0.050), dotted line 
(AE = 0.100) and long dashed line (AE = 0.500). 

which determines the position of the breather center as a function of time. 
The breathers are mobile and propagate unidirectionally with (virtually) con- 
stant velocity along the lattice up to an amount of disorder AE ~ 0.045. 
Generally, the larger the width AE the lower becomes this velocity. Finally, 
for an overcritical amount of disorder the breather performs oscillatory mo- 
tions around its starting lattice position. The frequency of these oscillations 
diminishes gradually with more enlarged AE and simultaneously the excur- 
sions from the starting site get continuously stronger confined. Despite the 
non-directed character of its motion even for a moderate amount of disorder 
AE = 0.05 the breather is able to travel as far as twenty lattice sites away 
from its starting position. Eventually linear Anderson-type disorder effects 
overwhelm more and more the nonlinear activation mechanism for breather 
motion through the excitation of the nonlinear pinning mode. 

Random double helix structure 

Concerning the impact of structural disorder we considered random distri- 
butions of the base pair spacings r n — r e [—A, A] and the twist angles 
9 n — 9q E [—A, A] with different mean standard deviations A. 

For representative results we averaged over several realizations of structural 
disorder. The results for the temporal evolution of the electron breather center 
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Fig. 6. Structural disorder and random double helix structure. The temporal behav- 
ior of the position of the electron breather center for different degrees of randomness 
in the radii r n and angles 9 n . Assignment of the line types to the increasing mean 
standard deviations ranging up to 10% of the respective mean value: Full line (1%), 
dashed dotted line (5%) and dashed line (10%). 

are depicted in Fig. 6. We conclude that structural disorder does not signif- 
icantly affect the mobility of the breathers. For growing degree of disorder 
merely the propagation velocity gets reduced. Even in the case of fairly strong 
disorder A = 10% (dashed line) the breather is still able to propagate steadily. 

We studied also the polaron mobility under the impact of combined static 
diagonal and structural disorder and found that the strong pinning effect gen- 
erated by static diagonal disorder prevails over the more harmless and merely 
velocity-reducing influence of structural irregularities. Particularly the last re- 
sult is in compliance with the findings in [41] where it was reported that 
hopping conductivity in proteins and DNA is influenced first of all by the 
molecular conformations and not by their sequences. 



5 Floquet analysis 



It is instructive to consider also the linear limit case of the coupled nonlinear 
electron-vibration system even though the results obtained throughout this 
paper bear essentially nonlinear character. The study of the linear limit system 
itself and its Floquet analysis can be applied as a reference and as a tool for 
obtaining useful information for the interpretation of the Floquet analysis of 
the fully nonlinear system. 
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5. 1 The linear system 



The linear limit system is derived from the Eqns. (12)-(14) for a = k = 
giving 

irc n = —Cn+i - c n _i (27) 
r n = -r n (28) 

Qnn-l — — ^ 2 Qnn-1 (29) 

In the first equation the terms E® c n , which in a homogeneous system re- 
duce to E c n , have been eliminated through the gauge transformation c n — > 
c n exp(— i Eq t/r). In the linear limit the electronic and vibrational degrees of 
freedom are decoupled. 

To construct the electronic state we substitute in Eq. (27), c n (t) = <p n exp(—iEt/r) 
with time-independent <p n yielding the linear stationary discrete Schrodinger 
equation: E(f> n — — n +i~ <t>n-i ■ With the Bloch mode solutions n = exp(i qn) 
we get then the dispersion relation: E = —e iq — e~ iq = —2 cosg. The linear 
electronic energy spectrum runs from —2 to 2. The polaron in the full system 
displayed in Fig. 2 possesses an energy value E = —2.064 which is slightly 
below the lower band edge E = — 2 assigned to a wave number q = and in 
the linear regime the frequencies are in the range w p = E/r G [—7.75,7.75]. 
Conclusively, the time scale for the evolution of the angular (radial) variables 
&nn-i { r n) represented by independent harmonic oscillators with frequencies 
1 and Q = 0.0842, respectively, is two orders (one order) of magnitude larger 
than the one of the electron. 

The linear analysis is also useful to understand the origin of many of the nu- 
merically calculated Floquet arguments for the fully nonlinear system. The 
values of the variables r n and 9 nn -i after a period T = 2n/w p = 2ttt/E are: 
r n (T) = exp(i2n /w p ) r n (0) and 6> nn _i(T) = exp(i 2 tt VL/w p ) 6> nn _i(0). There- 
fore, their corresponding Floquet arguments are given by 2 ir/w p and 2 7iQ/w p , 
for r n and # nn _i, respectively, which are attributed to angular positions in the 
complex plane with values around ±45° and ±3.9° for polarons with energies 
slightly below E = — 2 and with increased \E\ these angular positions shift to 
lower values on the unit circle. 



5.2 Results of the Floquet analysis 



The numerical integration of the system of tangent equations over one period 
of the polaron's oscillation T = 2 irr/\E\ allows us to obtain the actual Floquet 
matrix. Its eigenvalues can be written as exp(i ■?/>„) with ip n being the Floquet 
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arguments. The system is linearly stable provided all the Floquet arguments 
are real. As noted earlier in this paper within the range of parameters studied 
in this paper linear stability is assured, but this is not the only information 
that can be gained from the Floquet analysis. 

In order to relate the positions of the numerous eigenvalues of the Floquet 
matrix in the complex plane to the various degrees of freedom contained in the 
system it is elucidating to study first an anti-coupling limit in the electronic 
subsystem. Introducing an auxiliary tunable coupling parameter < /i < 
1, which corresponds to the transfer integral in the scaled system, the anti- 
coupling limit is obtained with \i = and the actually studied coupled system 
with ii = l. 

For the decoupled electron system (/i = 0) the value of the energy E can be 
easily derived with the help of the Eqs. (15) and (19), giving EQ n = — A; 2 |$| 2 $ n 
and r^°) = — /c|$ n | 2 , which leads to E — — 1 for the localized state $ n = 5 n Q 
and k = 1 as a typical value. Therefore, the r„-multipliers are in the vicinity 
of 7r/2, the # n „_i's lie at small angles at both sides of (1,0) and the c n 's are 
located at (1, 0) on the unit circle, in accordance with the analysis performed 
in 5.1. As we switch on the coupling ji and increase it up to the actual value 
H — 1, we can follow the evolution of the multipliers. Eventually, for fi = 1 
the r ra 's and # nri _i's have moved to positions with about half their initial 
angle values while spreading a little on the unit circle (the nn _i's to lesser 
extent though). The c n 's separate from (1,0) and invade the whole circle, as 
a consequence of the nonlinear coupling between the electron amplitudes and 
the vibrations. The exceptions are the pair at (1,0), corresponding to the 
phase mode, and another isolated pair close to this point. The latter pair is 
associated with the spatially antisymmetric pinning mode, responsible for the 
mobility of the polaron. Note that the linear approximation described in 5.1 is 
an approximation for the radial and angular variables. The Floquet multipliers 
for the decoupled and full systems, respectively are shown in Fig. 7. 

In order to illuminate further the role played by the vibrational modes with 
respect to the initiation of mobility of the polaron we cast the Floquet matrix 
F in the form of a supermatrix formed by submatrices: 

/ F F F i ^ 

1 cc 1 cr ± ct 



F F Ft 

± r c A rr ± r t 



^F 6c F 6r Fee ) 



(30) 



and the indices of the submatrices refer to the origin of their entries. From the 
arguments of the eigenvalues of the submatrices an/or combinations of them, 
it is possible to deduce additional information about the spectral features 
attributed to the various degrees of freedom contained in the coupled nonlinear 
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Fig. 7. The Floquet multipliers in the complex plane. Left figure: The linear limit 
system (27)-(29) at the anti-coupling limit \i = 0. Right figure: The fully coupled 
nonlinear system. 

system. The matrix F cc appears in the Born-Oppenheimer approximation. To 
recover the pinning mode of the complete system it suffices to consider only the 
subsystem consisting of the electronic part coupled to the radial motions for 
which the Floquet matrix reduces to the submatrix F C r = (F cc F cr ; F rc F rr ), 
which means that the mobility of the polaron is mainly a consequence of the 
coupling between the electron amplitudes and the radial variables and that 
the electron will not move within the Born-Oppenheimer approximation. The 
Figs. 8 and 9 

illustrate these statements. The left (right) part of the Fig. 8 shows the mul- 
tipliers of F rr {Fqq). In comparison with the linear case studied in the preced- 
ing section we note that while the radial components of the multipliers shifted 
their positions to lower angles and exhibit some spread the angular multipliers 
hardly alter their position. Hence, perturbations of the angular twist motions 
in the context of the coupled system resemble the behavior of the linear limit 
system. The left part of Fig. 9 depicts the eigenvalues of F cr in the vicinity of 
(1, 0). The pair of pinning mode eigenvalues is recognizable. In the right part 
of the figure the normalized pinning mode is portrayed. 

Moreover, it should be pointed out that for larger elongations of 9 nn -i, beyond 
the range where the linear expansion of <i„„_i is valid, there exists also an 
angular pinning mode as a consequence of the coupling between the electronic 
amplitudes and the angular variables. 

The fact that the disorder eventually prevents the polaron from being moved 
can also be explained with the help of the Floquet multipliers. When the 
degree of disorder is increased the pinning mode eigenvalues move progressively 
away from (1, 0) towards larger arguments ip and approach the rest of the F r 
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Fig. 8. Multipliers of F rr (left figure) and Fgg (right figure). 




Real part Sites 



Fig. 9. Multipliers of F cr in the vicinity of (1, 0). The pinning mode eigenvalues are 
recognizable. The right part of the figure shows the normalized pinning mode. 

eigenvalues. Fig. 10 displays the growth of the arguments of the pinning mode 
as a function of the disorder strength AE. Eventually the system is too far 
away from the instability due to the marginal mode, which appears when the 
pair of pinning mode eigenvalues collapses at (1, 0). This is easy to understand, 
as the disorder brings about a decoupling of the variables, favoring localized, 
stationary, isolated modes. 



6 Summary 

In the present work we have studied an electron transfer mechanism in DNA 
relying on the coupling of the charge carrying unit to the vibrational modes of 
DNA responsible for the formation of polarons and electron-vibron breathers. 
The steric structure of the double helix of A— DNA has been described in the 
context of the base-pair picture considering angular and radial vibrational 
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Fig. 10. Arguments of the Floquet multipliers of the low-frequency modes as a 
function of diagonal disorder AE demonstrating the growth of the argument of the 
pinning mode with increased amount of disorder so that the pinning mode eventually 
merges in the rest of the modes. 

motions of the base pairs. We have focused interest on the initiation of long- 
range and stable polaron and breather motion along the DNA structure when 
static and/or structural disorder inherent to any real DNA molecule is taken 
into account. We remark that the proposed polaron-like transport mechanism 
applies to electron hole conduction too. In the first part of the paper we have 
constructed exact stationary localized solutions (polarons and static electron- 
vibron breathers) with the help of a nonlinear map approach. We have found 
that the interplay of disorder and nonlinearity amplifies the degree of local- 
ization in comparison with cases when either disorder or nonlinearity alone 
is present. With concern to mobile electron- vibron breather solutions, whose 
motion is activated through a kick mechanism utilizing the pinning mode, 
with the current work we have demonstrated that they are fairly robust when 
the electron-vibration system is subjected to randomness. Generally, even for 
a moderate amount of disorder the breathers retain their localized shapes 
and support directed, long-ranged coherent electron transfer along strands of 
the bent double helix. In particular the last result suggests that DNA seems 
suitable for the design of (one-dimensional) functional nanostructures as ingre- 
dients of molecular nanoelectronical devices. Furthermore, we have observed 
that the electron breather is accompanied by vibrational breathers in the ra- 
dial as well as angular components, respectively causing local and temporal 
deformations of the traversed region of the double helix. Conclusively, efficient 
charge transport in DNA proceeds on the condition that the double helix un- 
dergoes structural changes which is exemplary for the interplay of structure 
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and function in flexible and adaptive biomolecules. 

Acknowledgments 

One of the authors (D.H.) acknowledges support by the Deutsche Forschungs- 
gemeinschaft via a Heisenberg fellowship (He 3049/1-1). J.F.R.A. acknowl- 
edges D.H. and the Institut fur Theoretische Physik for their warm hospitality. 
J. A. is grateful for a scholarship supplied by the Studienstiftung des deutschen 
Volkes. The authors are also grateful to the support under the LOCNET EU 
network HPRN-CT-1999-00163. 

References 

[I] D.D. Eley and D.I. Spivey, Trans. Faraday Soc. 58, 411 (1962). 

[2] C.A. Mirkin, R.L. Letsinger R.C. Mucic and J.J. Storhoff, Nature 382, 607 

(1996) ; A. A. Voityuk, N. Rosch, M. Bixon, and J. Jortner, J. Phys. Chem. 104, 
9740 (2000); F.C. Grozema, Y.A. Berlin and L.D.A. Siebbels, J. Am. Chem. 
Soc. 122 10903 (2000); Y.A. Berlin, A.L. Burin and M.A. Ratner, J. Am. Chem. 
Soc. 123, 260 (2001). 

[3] M. Ratner, Nature 397, 480 (1999). 

[4] D.B. Hall, R.E. Holmkin and J.K. Barton, Nature 382, 731 (1996). 
[5] D.B. Hall and J.K. Barton, J. Am. Chem. Soc. 119, 5045 (1997). 

[6] M.R. Arkin, E.D.A. Stemp, S.C. Pulver and J.K. Barton, Chem. Biol. 4, 369 

(1997) . 

[7] Y. Okahata, T. Kobayashi, K. Tanaka, and M. Shimomura, J. Am. Chem. Soc. 
120, 6165 (1998). 

[8] E. Meggers, M.E. Michel-Beyerle and B. Giese, J. Am. Chem. Soc. 120, 12950 

(1998) . 

[9] H.-W. Fink and C. Schonenberger, Nature 398, 407 (1999). 

[10] P. Tran, B. Alavi and G. Gruner, Phys. Rev. Lett. 85, 1564 (2000). 

[II] E. Braun, Y. Eichen, U. Sivan and G. Ben-Yoseph, Nature 391, 775 (1998). 

[12] D. Porath, A. Bezryadin, S. de Vries and C. Dekker, Nature 403, 635 (2000). 

[13] P.J. de Pablo, F. Moreno-Herrero, J. Colchero, J. Gomez Herrero, P. Herrero, 
A.M. Baro, P. Ordejon, J.M. Soler and E. Artacho, Phys. Rev. Lett. 85, 4992 
(2000). 

[14] J. Jortner, M. Bixon, T. Langenbacher and M.E. Michel-Beyerle, Proc. Nat. 
Acad. Sci. USA 95, 12759 (1998). 



22 



[15] Y.-J. Ye, R.-S. Chen, A. Martinez, P. Otto and J. Ladik, Sol. Stat. Comm. 112, 
139 (1999). 

[16] R. Bruinsma, G. Griiner, M.R. D'Orsogna and J. Rudnick, Phys. Rev. Lett. 
85, 4393 (2000). 

[17] Z.G. Yu and Xueyu Song, Phys. Rev. Lett. 86, 6018 (2001). 

[18] M. Hjort and S. Stafstrom, Phys. Rev. Lett. 87, 228101-1 (2001). 

[19] D. Ly et al, J. Am. Chem. Soc. 118, 8747 (1996). 

[20] J. Jortner, Proc. Nat. Acad. Sci. USA, 95, 12759 (1998). 

[21] E. Conwell and S.V. Rakhmanova, Proc. Nat. Acad. Sci. USA, 97, 4556 (2000). 

[22] D. Ly, L. Sanii and G.B. Schuster, J. Am. Chem. Soc, 121, 9400 (1999). 

[23] Z. Hermon, S. Caspi and E. Ben-Jacob, Europhys. Lett. 43, 482 (1998). 

[24] M. Peyrard and A.R. Bishop, Phys. Rev. Lett. 62, 2755 (1989). 

[25] L.V. Yakushevich, Quart. Rev. Biophys. 26, 201 (1993). 

[26] M. Barbi, S. Cocco and M. Peyrard, Phys. Lett. A 253, 358 (1999). 

[27] L. Stryer Biochemistry, (Freeman, New York, 1995). 

[28] S. Cocco and R. Monasson, J. Chem. Phys. 112, 10017 (2000). 

[29] T.D. Holstein, Ann. Phys. NY 8, 325, 343 (1959). 

[30] P.W. Anderson, Phys. Rev. 109, 1492 (1958). 

[31] N.F. Mott, Adv. Phys. 10, 107 (1961). 

[32] D. Bambusi and A. Giorgilli, J. Stat. Phys. 71, 569 (1993). 

[33] D. Hennig, J. Math. Phys. 39, 3568 (1998). 

[34] R.S. MacKay and S. Aubry, Nonlinearity 7, 1623 (1994). 

[35] D. Hennig, in preparation. 

[36] G. Kalosakas, S. Aubry and G.P. Tsironis, Phys. Rev. B 58, 3094 (1998). 

[37] N.K. Voulgarakis and G.P. Tsironis, Phys. Rev. B 63, 14302 (2001). 

[38] M. Peyrard and M.D. Kruskal, Physica D 14, 88 (1984). 

[39] S. Flach and C.R. Willis, Phys. Rep. 295, 181 (1998). 

[40] D. Chen, S. Aubry, and G.P. Tsironis, Phys. Rev. Lett. 77, 4776 (1996). 

[41] Y.-J. Ye and J. Ladik, Int. J. Quant. Chem. 52, 491 (1994). 



23 



